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We derive effective Hamiltonians for the fractional quantum Hall effect in n = and n = 1 
Landau levels that account perturbatively for Landau level mixing by electron-electron interactions. 
To second order in the ratio of electron-electron interaction to cyclotron energy, Landau level mix- 
ing is accounted for by constructing effective interaction Hamiltonians that include two-body and 
three-body contributions characterized by Haldane pseudopotentials. Our study builds upon pre- 
vious treatments, using as a stepping stone the observation that the effective Hamiltonian is fully 
determined by the few-body problem with N — 2 and A'^ = 3 electrons in the partially filled Landau 
level. For the n = case we use a first quantization approach to provide a compact and transparent 
derivation of the effective Hamiltonian which captures a class of virtual processes omitted in earlier 
derivations of Landau-level-mixing corrected Haldane pseudopotentials. 
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I. INTRODUCTION 

In a two-dimensional electron gas external magnetic 
fields form macroscopically degenerate groups of single- 
particle kinetic energy eigenstates known as Landau lev- 
els (LLs). Many-electron state degeneracies grow expo- 
nentially with system size when these Landau levels are 
fractionally occupied and electron-electron interactions 
and disorder are ignored. In this limit the zero tempera- 
ture chemical potential /i is constant over integer width 
intervals of the Landau level filling factor v, jumping be- 
tween single-particle eigenenergies at integer values of v. 
{f = N /Nll where N is the number of electrons in the 
gas and Nll is the number of single-particle states in 
each Landau level.) The fractional quantum Hall effeclP 
is a consequence of jumps in /i at non-integer values of 
V, and therefore can occur only in interacting electron 
systems. 

Because of the exponential degeneracy of the many- 
body ground state when interactions are neglected, it is 
not possible to understand the fractional quantum Hall 
effect by treating interactions as a weak perturbation. 
Instead, the problem of interactions in systems with par- 
tially filled Landau levels has traditionally been simpli- 
fied by allowing occupation numbers to fluctuate only 
within the partially filled level. This projection of the 
interaction Hamiltonian onto a single Landau has a long 
history in theories of two-dimensional electron systems, 
and was first employecP^ even prior to the fractional 
quantum Hall effect's discovery.^ It is strictly justified 
as a low energy theory, only when the interaction energy 
per particle is small compared to the energetic separa- 
tion between the Landau levels. In this article we derive 
effective Hamiltonians which account for corrections to 
the projected Hamiltonian that are valid to leading or- 
der in interaction strength. These corrections account for 
quantum fluctuations in otherwise empty and full Landau 
levels, and are therefore normally referred to as Landau 
level mixing corrections. 

Because of its non-perturbative character, the problem 



of interactions in a system with a partially filled Landau 
level has been a rich source of unique correlated electron 
phenomena, including fractional and non-Abelian quasi- 
particle statistics,^ and electron-hole pair superfiuidity.'^ 
We limit our attention in this paper to the case of a two 
dimensional electron system with parabolic bands. The 
semiconductor quantum well systems in which the frac- 
tional quantum Hall effect has most often been observed 
are well described by such a model. The small parameter 
on which our analysis is based is the ratio of characteristic 
interaction and kinetic energy parameters in the strong 
magnetic field limit of a two-dimensional parabolic band: 



where ujc = eB/m*c is the cyclotron frequency, Iq = 
\fhcJeB is the magnetic length, m* is the parabolic 
band effective mass, and e is the low-frequency dielec- 
tric constant of the environment hosting the two dimen- 
sional electron system. (Hereafter /i = 1.) Note that k, 
varies as I/^/B^ reaching its smallest values at the largest 
fields. The fact that the projected Hamiltonian is able 
to provide an adequate description of most properties 
of systems with fractionally filled Landau levels is, at 
first sight, somewhat surprising. In electron-doped gal- 
lium arsenide (GaAs), for example, k ~ 0.4 even at the 
highest achievable steady magnetic fields, B ^ 40T. In 
hole-doped GaAs,'^ AlAs,'^ and in the recently studied 
ZnO heterostructure^ effective masses are larger, reduc- 
ing the cyclotron energy and increasing k further, as sum- 
marized in Table |T] Landau level mixing is also expected 
to be substantial in graphen^ and in silicon quantum 
wells. Additionally, some of the most interesting frac- 
tional quantum Hall states occur in higher levels, and 
therefore at weaker magnetic fields and hence larger k 
values. Our goals in this paper are to shed light on why 
the influence of LL mixing on the fractional quantum 
Hall effect is often modest, and to make progress in un- 
derstanding its role when it is essential. 
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There are many specific motivations for the study of in- 
teraction induced LL mixing. One is to attempt to bring 
theory and experiment into closer quantitative agreement 
in cases where the quahtative picture is ah'eady under- 
stood. It is widely recognized, for example, that theo- 
retically predicted FQHE gaps are invariably larger than 
the experimentally measured ones.'^^' Although part of 
the discrepancy can be attributed to disorder and to fi- 
nite quantum well widths L L mixing is also expected 
to play an important role.'^^ISHni 

Another motivation is to study the emergence of phases 
which would otherwise be unstable. Of particular inter- 
est is the potential role of LL mixing on the stabilization 
of the Moore-Read Pfafhan statej^ generally believed to 
describe the incompressible state observed at filling frac- 
tion V = 5/2.^^-^° Indeed the Pfaffian and its particle- 
hole conjugate, the anti-Pfaffian, are distinct phases,'2il22l 
which are energetically degenerate in the absence of LL 
mixing. Which of these two states appears experimen- 
tally is completely determined by the particle-hole sym- 
metry breaking terms that Landau level mixing gener- 
ates. This state selection property for half-filled LLs ap- 
plies to any incompressible state that is not particle-hole 
invariant. Many numerical studies support the view that 
a state of the Moore-Read type is favored atv — ^/2^^^ 
but the detailed form of the LL mixing would determine 
which of the two particle-hole mirror states is preferred. 
It is still unsettled which one of these two states is se- 
lected by the LL mixing induced by pure Coulomb in- 
teractions. One study has found the Pfaffian state to 
be favored,'^^ while another has found the anti-Pfaffian 
to be favored,!^ although their detailed account for the 
LL mixing was different. Depending on the details of 
the confinement, the intersubband LL mixing could also 
play an important role,065 specially in wide wells, where 
experimental studies have highli ghte d its influence of on 
the stability of the v = 5/2 statel^ 

Additionally, differences in gap sizes between filling 
factor n+i> and n+l—h' (or 2n+iy and 2n+2—iy when spin 
is an active degree-of- freedom), observed conspicuously 
in GaA^2l! and in large k systems like AIAJ^SS and hole- 
doped GaAs,!^ require particle-hole symmetry breaking 
and hence LL mixing. Finally, particle-hole asymmetries 
in critical densities for the emergence of charge density 
wave states at low particle or hole densitie^^ also reflect 
LL mixing. 

In this article we will construct an effective Hamilto- 



TABLE I: Relative strength of Coulomb and cyclotron 
energies in GaAs, AlAs^ and ZnO^ {B is measured in 
Tesla). 





m* /mo 
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electron GaAs 


0.069 


13 


2.6/v 




hole GaAs 


0.39 


13 


14.6/^ 


/B 


electron ZnO 


0.29 


8.5 


16.7/^ 


/B 


electron AIAs 


0.46 


10 


22.5/^ 


/B 



nian which accounts for LL mixing to leading pertur- 
bative order in k. This many-body Hamiltonian when 
solved exactly will be able to predict energies, in units of 
Wcj correctly to order k^, and the projected many-body 
wavefunctions into the Landau level of interest to order 
K. It is interesting to note that the new energy scale ap- 
pearing in our analysis, WcK^, is independent of magnetic 
field strength (unlike the dominant interaction scale of 
the FQHE, uicH, which grows as VB), and it is twice the 
effective Rydberg of the parabolic band system: 



UlcK 



(2) 



One of the earliest studies to account numerically for 
LL mixing in the second Landau level was performed by 
Rezayi and Haldane.^^ Analytical perturbative studies 
of LL mixing due to interactions were first carried out 
in the lowest Landau level by Murthy and Shankar.^i^ 
These studies were subsequently extended to the sec- 
ond Landau level by Bishara and Nayakji^Sl employing an 
analysis similar to the renormalization group (RG) for 
fermions. The present work complements these earlier 
papers by presenting new methods of derivation, adding 
some new results, and correcting some previous results. 
Our work has been developed essentially in parallel with 
two recent studies by Peterson and NayakjISI and Simon 
and Rezayi, ^'^ and our findings are largely in agreement 
with these two studies. A key observation in our ap- 
proach is that the effective many body Hamiltonian can 
be constructed by solving N — 2 and A^ = 3 few-body 
problems, which we use to simplify some derivations. 
We compute two-bodj'^^'^ and three-bodjE^ElIlZl gen- 
eralized Haldane pseudopotentials for these interactions, 
which can be incorporated into many-body numerical di- 
agonalization studies. Although the values we list here 
for these pseudopotentials are specialized to the case of 
2D Coulomb interactions, we have derived analytic and 
semianalytic formulae for all the pseudpotentials valid 
for any rotationally invariant interaction. Pseudopoten- 
tial parameter values for more realistic interaction mod- 
els which account for finite quantum well widths can be 
conveniently computed from these expressions. 

Our paper is organized as follows. Section [H] provides 
a compact derivation of the effective Hamiltonian in the 
lowest Landau level {n — Q LL) that is based on a first 
quantization formalism, and is valid for both bosons and 
fermions. In Sec. IIIII we construct the effective Hamilto- 
nian for a partially filled first excited Landau level (n = 1 
LL), and compute its two-body and three-body general- 
ized Haldane pseudopotentials. In Sec. IV we derive an 



effective Hamiltonian valid for 1 < < 2 that is valid for 
the special case of maximally polarized electronic states 
in which the majority spin state is full and the minor- 
ity spin state is partially occupied. We also compute 
the two and three-body generalized Haldane pseudopo- 
tentials appropriate for these effective Halmitonians. In 
Sec. [V] we summarize our findings and present conclu- 



3 



sions. We have relegated discussions of some effective 
interaction properties and calculation details to a series 
of appendices. 



II. LOWEST LANDAU LEVEL 

A. Many body Hamiltonian to order 

We consider a two dimensional electron system sub- 
jected to a perpendicular magnetic field B = —Be^. 
Measuring all energies in units of the cyclotron frequency 
and all lengths in units of the magnetic length, the single- 
particle spectrum of a non-interacting disorder-free spin- 
less system consists of discrete Landau levels with ener- 
gies. En = n + 1/2.'^^'^ The non-interacting Hamiltonian 
(including the Zeeman energy contribution) and the full 
Hamiltonian including interactions are given respectively 
by, 



1/2 - gcrt 



and 



(3) 



(4) 



where fii + 1/2 is the kinetic energy operator of particle 
i, g = gsTTi* /2me with me the mass of the electron in 
vacuum and gs the effective g-factor of the host mate- 
rial, and Vij = l/\ri — rj\ is the dimensionless Coulomb 
potential. 

Unless the filling fraction v is an integer, the non- 
interacting many-body eigenstates are degenerate. From 
degenerate state perturbation theory, the eigen-energies 
can be determined to order k by projecting the interac- 
tion term onto the degenerate manifold of non-interacting 
eigenstates with energy Eq, 



Hi=Eo 



i<j 



(5) 



"Hi is the Hamiltonian commonly employed to study the 
fractional quantum Hall effect. Employing conventional 
degenerate second order perturbation theory, the correc- 
tion to the next order, in energies and k in the pro- 
jected wavef unctions, can be obtained from the effective 
Hamiltonian, 



'H2 — Hi 



k<l 



-V^VkiV, (6) 



where V is the projector into the degenerate non- 
interacting ground state manifold of energy Eq, and 
V± = 1 — V, is the projector onto its orthogonal com- 
plement. In the absence of Zeeman energy, ■H2 is the 



first quantization version of the Hamiltonian considered 
in Ref. 17, 

In the case of the n = LL, the degenerate mani- 
fold would be a subspace of the lowest kinetic energy 
eigenspace with a definite projection of the total spin 
along the z-axis, S^, and Zeeman energy Eq ~ —2gSz. 
Nevertheless, the energy denominator in Eq. Q only in- 
cludes the kinetic energy difference between the virtual 
excited states and the states in the degenerate manifold. 
This is a consequence of the conservation of S'^, because 
the virtual excitations produced by the interactions Vij 
and Vki do not change it, and therefore, the Zeeman en- 
ergy disappears from this energy denominator. Conse- 
quently, for any 5*2 , we can write the second order correc- 
tion to the effective Hamiltonian when the lowest-energy 
n = LL is partially filled, as. 



— T~Li 



k<l 



n 



(7) 



with fi = f^i- There are three possibilities for the two 
pairs of indices (i,j) and (fc, /) appearing in the sum in 
Eq. ([7]); they can share both particle indices (i.e. i = k 
and j — I), they can share only one particle index while 
the other two are distinct, or they can share no indices. 
The last possibility does not contribute to because 
the projection operators force each virtually excited par- 
ticle to decay back into the lowest LL after both interac- 
tions act. The other two possibilities are non- vanishing 
and yield respe ctively two- and three-body effective in- 
teractions .E^HIl Below we address the two-body effective 
interactions first. 



B. Two-body interactions 

The effective interactions implied by Eq. ([t]) are in- 
dependent of the state of the many body system in the 
n = LL. We can therefore, without loss of generality, 
determine the effective Hamiltonian by considering only 
the few body N — 2 and = 3 cases. The two body 
interactions can be written as. 



V^'' = KVV12V - K^Vvi2V_L-V_LVi2V. (8) 

n 

We have explicitly verified that this interaction is identi- 
cal to the the two body interaction in Ref. 17,, In particu- 
lar, after translating V^'' from Eq. ([8| into its equivalent 
second quantized version, the piece for which both parti- 
cles are virtually excited into higher Landau levels corre- 
sponds to the interaction SHqq in Eq. (26) of Ref.lTTl and 
the piece of Eq. ([8 ) for which only one particle is excited 
into a higher Landau level corresponds, after normal or- 
dering is performed, to the two body part obtained from 



SH^Q in Eq. (30) of Ref. HH 
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The translational and rotational invariance of Eq. ([8]) 
permits its decomposition into Haldane pseudopotentials. 
The virtual excitations are most easily analyzed by de- 
composing the degrees of freedom into center of mass, R, 
and relative, r, coordinates. The total kinetic energy is 
then the sum of relative and center of mass kinetic ener- 
gies, h = hi+h2 = hr + fifj. Because the interaction acts 
only on the relative coordinate, only fir enters the energy 
denominator. We write the two-body Haldane pseudopo- 
tentials as Vj^ = Kyi(m) -|-K^F2(m ) where Vi (m) are the 
standard Haldane pseudopotentiala^^I^^I^^I^ given by 



Vi{m) = r(0,m|i;i2|0,m)r 



r(m+l/2) 
2m! 



(9) 



Here |n, m)r is an state with the relative guiding center 
and kinetic energy quantum numbers m and n respec- 
tively. Note that the form of the relative motion eigen- 
states differ from the familiar single-particle eigenstates 
only because of the difference between relative motion 
and single-particle motion magnetic lengths £r = V^£o- 
For V2(m) we obtain: 



{n, m + n|tii2|0, m)r'\' 



4(m + 1) 



1 1 3 3 

^' ^' 2 ' 2 . 2 

2,2,m-|-2' 



(10) 



with 4i^3 the generalized hypergeometric function. Val- 
ues of these pseudopotentials are presented in Table [Hj 
together with the usual Haldane pseudopotentials for 
comparison. These pseudopotentials agree with those 
computed in Refs.l^andfi5l Additionally very closely re- 
lated numbers have been computed previously in a strong 
magnetic-field expansion of the spectrum of hydrogenic 
atoms.l^ V2{m) equals 7rao,m/4 where ao,m values are 
specified by Eq. (15) and Table II of Ref 

Our results for the leading order Haldane pseudopo- 
tential corrections differ from those presented in Table I 
of Ref. im The reason for this discrepancy is that the set 
of virtual processes in which only one particle is excited 
into higher Landau levels was omitted in Ref. HH This 
discrepancy has been recently solved in Ref. [9l It was 
shown there that, in the RG language, these additional 
processes arise from keeping track of the correct normal 
ordering of the three-body interactions which give rise to 
additional contributions to the two-body interactions. 

Since we have made no explicit reference to the statis- 
tics of the particles involved, our results apply equally 
well to fermions and bosons. The difference between 
the two cases is only in the constraint imposed by quan- 
tum statistics on the allowed states, which implies that 
odd (even) m pseudopotentials are associated with spin 
triplet wavefunction for fermions (bosons) and spin sin- 
glets for bosons (fermions). The case of a partially filled 
n > LL, that we will discuss later in Sec. |III[ is rele- 
vant only for fermions since partially filled higher Landau 



level states are not non-interacting ground states in the 
bosonic case. 



C. Three-body interactions 

In order to derive the full Hamiltonian including the 
three body terms, it is sufficient to consider a = 3 few 
particle problem. Three-body interactions follow from 
the terms in Eq. ([7]) where only one index is shared be- 
tween the pairs («, jj and (fc, I) while the other two indices 
are distinct, e.g. i — k but j ^ I. This contribution can 
be written as 



3b 



-K^y2nsVvi3V±lv±vi2vns, (ii) 



seS3 



where the sum is over the six permutations of three ob- 
jects and Us is the associtated unitary permutation oper- 
ator. Since bosonic and fermionic states are both eigen- 
states of the permutation operator Hgl^*) = ilvE"), it fol- 
lows that 



($|V3''|^') = -6k2($|7'wi3'P_l--P_lwi2-P|*), (12) 

n 

for arbitrary states {j^E'), |<i>)} with the same parity un- 
der permutations. Generalized Haldane pseudopotentials 
for A'^-body interactions have been thoroughly discussed 
for spinless particles by Simon, Rezayi, and Cooper in 
Ref. 133 and for spinful particles by Davenport and Si- 
mon in Ref. |47l We specialize hereafter in the case of 
fermions. There are two possible values for the total spin 
of three particles, namely S = 3/2 and S = 1/2. The 
spatial wavefunction for S = 3/2 must be fully antisym- 
metric, while the one for 5=1/2 has mixed symmetry. 
For S = 3/2 we employ the fully antisymmetric wave- 
functions for three particles in the lowest LL constructed 
by Laughlin in Ref. |5TJ whose polynomial part is 



kl ~^ ^3/2 



(zg + izfc)^' - {zg - izbf^ 
2i 



(13) 

with Z^f = 23'+2fc+i[7r3(3/ -I- fc)!fc!]i/2, k>0,l>l, and 



2 f Zi+ Z2 

3 



Z3 



Zb = 



Zl - Z2 



(14) 



These wavefunctions form a complete orthonormal ba- 
sis for the relative internal states of three particles, 
i^k'H'^lfi') = 5k',k5i',i. The relative orbital angular 
momentum of the states is m\^]!j') = {2k + 3^)|\E'^;^). 
Therefore, it follows from the rotational invariance of in- 
teractions that the generalized Haldane pseudopotential 
matrix satisfies, 
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TABLE II: Conventional two-body Haldane pseudopotentials in the n = LL (Vi from Eq. (|9| 

of their leading pertubative corrections (V2 from Eq.(lO)). 



and the coefficients 



m 





1 


2 


3 


4 


5 


6 


7 


8 


9 




0.8862 


0.4431 


0.3323 


0.2769 


0.2423 


0.2181 


0.1999 


0.1856 


0.1740 


0.1644 


V2 


-0.3457 


-0.0328 


-0.0112 


-0.0055 


-0.0033 


-0.0022 


-0.0015 


-0.0012 


-0.0009 


-0.0007 



TABLE III: Coefficients of the leading S = 3/2 (V3/2 from Eq. ([Ts])) and S = 1/2 {V1/2 from Eq. ^) three-body 
Haldane pseudopotentials in lowest LL. For S = 3/2 there are two states with total angular momentum m = 9 so 
the pseudopotential is a 2 x 2 matrix. For = 1/2 there are two states for 4 < m < 6 and Haldane pseudopotentials 
are matrices, in these cases the listed pairs have been orthonormalized by rotating only the state with cr = 2, as 

described in Appendix [A| 



m 


3 


5 


6 


7 


8 


9 


(fc,0 

V;/2{k'l',kl) 


(0,1) 
-0.0181 


(1,1) 
0.0033 


(0,2) 
-0.0107 


(2,1) 
0.0059 


(1,2) 
-0.0048 


(0,3) (3,1) 

/ -0.0049 -0.0007-1 
I -0.0007 0.0052 ) 


m 


1 


2 


3 


4 


5 


6 


{a,k,l) 
Vi/2{a'k'l',akl) 


(1,0,0) 
-0.0345 


(2,0,0) 
-0.0540 


(1,1,0) 
0.0425 


(1,0,1) (2,1,0) 

/ -0.0343 -0.0025 \ 
l -0.0025 0.0075 ) 


(1,2,0) (2,0,1) 

/ 0.0277 0.0067 \ 
I 0.0067 -0.0176 j 


(1,1,1) (2,2,0) 

/ -0.0119 -0.0050-1 
I -0.0050 0.0102 J 



(K'H^^'IK?) = '^'S2k'+3l',2k+3lVy2ik'l',kl). (15) 

In this way we find the S — 3/2 pseudopotentials listed 
in Table [Till These pseudopotentials agree with those de- 
rived in Refs.[^and.42, They are also in close agreement 
with the pseudopotential differences obtained in Ref. [35l 
Further details on the properties of these states and the 
derivation of these pseudopotential values can be found 
in Appendix \K\ 

To construct the 5=1/2 states we follow the approach 
of Davenport and Simon.l^S We start with two primitive 
polynomials which are antisymmetric only under permu- 
tations of variables 1 and 2, and thus have the symmetry 
of the Young tableaiP^l 



1 


3 







(16) 



for the fully symmetric polynomials a different one from 
that of Ref. ,47. Our basis for the symmetric polynomi- 
las is the bosonic analog of the fermionic wavefunctions 
of Ref. 1511 By multiplying this bosonic wavefunction by 
the primitive polynomials, /3cr, we obtain the polynomial 
part of the = 1/2 fermionic wavefunctions. 



yl/2 



(4 + 4)' 



{Za + iZbY^ + {Za - iZb) 



311 



with a = {1,2}, k > 
tion factor for ct = 1 



(18) 

0, / > 0. The normaliza- 
statcs is given by 



^fc/ V(l + <5i,o)(2fc + 2 + 30, and for a 



2 it is 



Zli'' ^(1 + Si,o)[ik + 1 + 3/)(fc + 2 + 30 + (fc + l)(fc + 2)]. 

The spin part of the 5* = 1/2 wavefunction has the 
symmetry of the Young tableau conjugate to (16), thus 
the fully antisymmetric wavefunction 



These primitive polynomials are, 



Pi = Zb, P2 = ZaZb. 



(17) 



The most general polynomials with the symmetry of this 
tableau are then obtained by multiplying the primitive 
polynomials by the most general fully symmetric trans- 
lationally invariant polynomial. We choose as a basis 



(19) 



where A = (1 — H23-I-Hi2H23)/a/3, is the partial antisym- 
metrization operator between f and \, particles, and Wij 
is the unitary operator corresponding to elementary per- 
mutations between i and j. The states are linearly 
indepenent and complete but not orthogonal: 



o4fe'+6/' + l 
if 

yl/2 71/2 



[k\{k'+il')\5v,i+i5k,k'+i{l + Sifi) + k'\{k + 30!<J/,i'+i4',fc+2(l + Svfl)]- 



(20) 
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The internal motion orbital angular momentum of these 
states is m\'^lZ) = + 2fc + 301*^^,)- Thus the gener- 
alized Haldane pseudopotentials satisfy 



(K'l'l'l^^'KZ) - ^'S,,+2k'+3l'.a+2k+3lV,/2{a'k'l', <jkl), 

(21) 

and have the leading values listed in Table |III[ These 
pseudopotentials are also in agreement with those de- 
rived in Refs. and SH where the two leading three-body 
pseudopotentials for S = 1/2 were computed. Further 
details on these states and the 5 = 1/2 three-body pseu- 
dopotentials can be found in Appendix \K\ 



III. n = 1 LANDAU LEVEL 
A. Many body Hamiltonian to order 



Landau level. These produce only a well known over- 
all constant shift of the single-particle energies, that is 
equivalent to a change in the chemical potential.^^ 

The non-vanishing two-body terms satisfy a kinetic 
energy balance condition, + uq =713-!- 714, and can 
take values + 714 = {2,1,0} in Eq. (23). We em- 



ploy the three possible values of this incoming kinetic 
energy as labels for the three allowed terms, labeled as 
{V2^(2), V2''(1),V2''(0)}, and find that: 



flu 



Wl2,65l'56,34 Z C{^C^^C3fj,C4u, 



1---6 



(24) 



The first quantization analysis presented in section |TT] 
for a partially filled lowest LL is cumbersome in the case 
of a partially filled n = 1 Landau level because of the 
need to account for virtual excitations of the full n = 
LL. In this section we therefore use a second quantization 
approach. By following a path entirely analogous to that 
of section II A we arrive at an effective Hamiltonian to 
order which acts in the n = 1 LL. 

The first order term is as usual simply the projection 
of the Hamiltonian onto the partially filled level, 

Hi = i?0 + - X Wi2,34 4miiy4m2M'^l'"3A'Cim^^, (22) 

where c]^^^ is a fermion creation operator in LL n with 
guiding center quantum number m, fj, runs over all the 
fermion flavors, and Eq is the non-interacting energy in- 
cluding the Zeeman term. For generality we assume /- 
flavors (i.e. for spin s, / = 2.s + 1). To order the 
effective Hamiltonian is 



where the indices of the creation/annihilation operators 
are understood to be on the n = 1 LL, 6(nc,^nQ) is a 
function that restricts ~\- > 3 and > l,ng > 1, 
and S{n5, ng) restricts — uq — 0. The term V^^{1) can 
be split into three terms. 



V2''(l) = Vf(l) + V2''(l) + Vf(l), 



K (1) = 2^ Wl5,64W26.53 ci^C^^C3^C4i,, 

= 2^ Wi5, 36^62,54 ci^ci^C3pC4i., 

Vf''(l) = X(wi5,46«62,35 + W51,46^'26,35) X • • • 



Ci^C2„C3^C4i,, 



(25) 



'H2 — l-il — X 1'l2,34l'56,78 



1---8 



1 



P±cl^cl^c-i„C4),Po, (23) 



fi- Ni 



where the integers abbreviate single-particle kinetic and 
guiding center quantum numbers (e.g. {1 ^ nimi}), A^i 
is the number of particles in the partially filled n = 1 
LL, Pq is the projector into the many-body eigenspace 
with fi — Ni, and P± is the projector into its orthogonal 
complement ^ 

The classification of all the interaction terms arising 
from this Hamiltonian is a lengthy bookeeping exercise 
that we describe in Appendix [B] There are no four body 
terms, and the one body terms, which we do not com- 
pute, account only for exchange interactions with the full 



where r^n^^ng) restricts — and ng > 1. Each of 
these two-body interactions has a simple diagrammatic 
representation. In particular V^^(2) and V^''(0) corre- 
spond to particle-particle ladder diagrams, and V^^(l) is 
a screening diagram, hence the flavor multiplicity factor. 
V^''(l) is a particle-hole ladder diagram, and V^''(l) is a 
vertex correction diagram as summarized in Fig. [ij To 
get a feeling for the signiflcant additional complication of 
the n = 1 LL effective interaction calculation, note that 
from all of these terms only the particle-particle ladder 
type diagram is present in the n — LL. When our first 
quantization formulation approach was used for the two 
body interactions in the n — LL case, the problem 
was reduced to the simplicity of a one-body second-order 
perturbation theory calculation. 

By comparison, the three-body interactions are sim- 
pler. They can be combined into a single term, 
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x(2) 










2fj, 


lu 


5u 




6/i 








3(1 


4z/ 



x(2/) 



5a 




6cr 



the two-body Haldane pseudopotentials, we consider the 
2(1 N = fNLL + 2-body problem in which the lowest Landau 

level is completely full with fNi^^ particles and only two 
particles are in the n = 1 LL. The state describing the 
orbital part of the relative motion of two particles in the 
'^^ n = 1 LL, |m)i , can be written as/^ 



x(4) 



x(4) 



lu 




2n 


lu 




2/i 


lu 




(1) du 




5// 


vf(i) 








6u 




6(1 






5u 


4u 






'lu 




3/< 


4i/ 





2/i 



3/i 





X 


(4) 


2fi 




lu 






7u 


5fi 




5u 



3A 



4A 



FIG. 1: Diagrams representing two- and three-body 
interactions at order in the n = 1 LL. The symbol 
X (fc) appearing above each diagram represents its 
multiplicity, i.e. there are k different contractions of the 



operators in Eq. ( 23 ) which give rise to the same 
diagram after relabeling dummy indices. The indices of 
the fermion lines are consistent with those appearing in 
Eqs. (psl) and (p6|. 



E 

1--7 



^'12, 57^73, 46 



1 

nj — 1 



-1 Cg^;^ C4 A C5^ C6iy , 



(26) 

where again the Landau level indices of the operators 
are understood to have n — I and uj{nj) restricts n-j ^ 
1. The three-body interaction has a simple diagramatic 
representation depicted in Fig.[l] There is a compact first 
quantized version of Eq. ( |26[ ), which for three particles 
in the n = \ LL is 



V' 



i3& 



(27) 



where the notation is the same as that of Eq. (11), ex- 
cept that V is the projector into the eigenspace with 
n = ni +n2 +?T.3 =3, and V±_ into its orthogonal comple- 
ment. This first quantized version significantly simplifies 
the evaluation of three-body Haldane pseudopotentials 
which we will discuss later on. 

V2''(2), Vf{l), Vl\l), Vf (1), V2^(0), and are 
separately rotationally and translationally invariant, as 
discussed in Appendix [Cj This property guarantees that 
the interactions can be represented by Haldane pseu- 
dopotentials as in the n = LL case. To compute 



\m)i = a\a\\Q, m)r, 



(28) 



where |m, 0)^ is the state of relative angular momentum 
m, of two particles in the n = LL discussed next to 
Eq. (|9]), and aj is the kinetic energy raising operator for 
particle i. Note that |m)i is an entagled state of rela- 
tive and center of mass quantum numbers as discussed 
in Appendix [Cj The u sual Haldane pseudopotentials for 
the n = 1 LL areP^^ 



Vi(m) = i(m|t)i2|7Ti)i = 

r(m + 1/2) (m - 3/8)(m - 11/8) 
2m! (m- l/2)(m- 3/2) ' 



(29) 



The Haldane pseudopotentials to order are listed in 
Table |IV[ Our pseudopotentials are in agreement with 
those obtained in Refs. [HI lU and SSI with small dis- 
crepancies pressumably arising from numerical error. We 
believe our pseudopotentials are essentially free from nu- 
merical errors because we have converted the effective in- 
teraction expressions to first quantization, as dicussed in 
Appendix [cj which allows for very efficient calculations. 
We have explicitly verified that our two body interac- 
tions are equivalent to those of Ref . JL More specifically, 
the sum of our ladder type interactions from Eq. (24 1, 
namely V^''(2) -I- V^''(0), is the same as the sum of the 
BCS interaction of Eq. (19), with the two-body interac- 
tion appearing in the last line of Eq. (26) in Ref. [91 And 



the sum of the interactions appearing in Eq. (25), namely 
Vf (1) + Vfc2''(l) + V^^l), is identical to the sum of the 
ZS and ZS' interactions appearing in Eqs. (17) and (18), 
with the remainding two-body interactions appearing in 
the second and third line of Eq. (26) in Ref. [HI 

It is interesting to note from Table |IV| that the lead- 
ing two-body pseudopotentials in the n — 1 LL, namely 
those with m < 3, are dominated by the screening in- 
teraction V^^(l). The contributions from the remaining 
interactions to these pseudopotentials nearly cancel. In 
Appendix [c] we discuss further properties of Vf''(l), and 
explicitly show that it is equivalent to the term in the 
RPA approximation for the statically screened potential 
in the presence of a completelly filled n — LL. 

The three-body states of the TV = fN^^+S body prob- 
lem in the n = 1 LL are mapped from those in the n — 
LL by raising the kinetic energy of the three particles, 



(30) 
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TABLE IV: Conventional two-body Haldane pseudopotentials in the n 
their leading perturbative corrections for a completely filled spin- 1/2 n = 



1 LL (Vi from Eq. 



29) and coefficients of 
^ 2), 



LL (from Eqs. ([241 and ^ with / 
which describe the physics for filling factors v in the interval 2 < < 4. The last row is the sum of all the 
perturbative corrections, i.e. V^'' = V^''{2) + Vf (1) + (1) + (1) + V^'>{0). 



m 





1 


2 


3 


4 


5 


6 


7 


8 


9 


Vi 


0.6093 


0.4154 


0.4500 


0.3150 


0.2635 


0.2322 


0.2101 


0.1935 


0.1803 


0.1696 


V'"'(2) 


-0.0903 


-0.0347 


-0.1235 


-0.0241 


-0.0110 


-0.0064 


-0.0042 


-0.0030 


-0.0022 


-0.0017 


vf(i) 


-0.3930 


-0.2038 


-0.1981 


-0.1119 


-0.0535 


-0.0235 


-0.0098 


-0.0039 


-0.0015 


-0.0006 


Vf(l) 


0.0247 


0.0706 


0.0803 


0.0186 


-0.0031 


0.0091 


-0.0004 


0.0029 


0.0007 


0.0011 


V?(l) 


0.0750 


-0.0475 


0.0870 


0.0141 


-0.0128 


-0.0167 


-0.0139 


-0.0104 


-0.0077 


-0.0058 










-0.0276 


-0.0023 


-0.0006 


-0.0003 


-0.0001 


-7 X 10~^ 


-5 X 10~^ 


-3 X 10" 


y2b 


-0.3836 


-0.2155 


-0.1818 


-0.1056 


-0.0810 


-0.0377 


-0.0285 


-0.0146 


-0.0108 


-0.0070 



With this construction the orthonormality of the S — 3/2 
states in n = LL and Eqs. ( 15 ),( 20 ), and (pTl) are imme- 



diately extended to the n = 1 LL three body states. The 
three body Haldane pseudopotentials we obtain are listed 
in Table |V] Further properties of these states and of the 
corresponding pseudopotentials are given in Appendix [X] 
The three-body Haldane pseudopotentials in the n = 1 
LL are in agreement with those computed in Refs. ,9 and 
S^l There is a discrepancy with those obtained numer- 
ically in Ref. |33] which is likely to arise from the errors 
associated with the finite size effects and Landau level 
trunctaion present in such study. 



IV. PARTIALLY FILLED SPIN POLARIZED LLL 

Another instance for which a useful effective Hamil- 
tonian can be obtained using the line of reasoning pre- 
sented in the previous sections, is the case of a Landau 
level for which spin is a good quantum number, and fo 
spin states are completely filled {fo < f) while the re- 
maining f — fo spin states are partially empty. This ap- 
proach is useful for electrons, for example, in addressing 
those states at filling factors in the interval 1 < v < 2 
for which all f = states are occupied, i.e. for max- 
imally spin-polarized states. Although we could use the 
effective Hamiltonian discussed in Sec. [H] in this filling 
factor range, it is useful to derive an effective Hamilto- 
nian which acts only on J, degrees of freedom, since the 
conservation of Sz prevents the f spins from participat- 
ing in the low energy dynamics. We emphasize that the 
ground state in the filling factor range 1 < < 2 is not al- 
ways maximally spin-polarized, so this approach cannot 
always be used to describe the ground state. 

We can construct an effective Hamiltonian for the par- 
tially empty flavors in which the interaction is exactly the 
same as the one discussed in Sec. |TT] except that there is 
an additional two-body interaction of the screening type, 
in which occupied flavor electrons are virtually excited to 
higher Landau levels in a completely analogous manner 
to V^^{1) from Eq. (25), with an analogous diagramatic 
representation as that appearing in Fig.[l| except that its 



multiplicity will be given only by the /o occupied flavors. 
This additional contribution to the effective interaction 



is: 



Wl5,64t'26,53 C^j,C2^,C3^C4y , 



(31) 

where the orbital indices of the operators are understood 
to be in the n = LL, and the spin indices run over the 
partially empty f — fo flavors only, i.e. all electrons have 
minority spins when only spin provides a flavor label. 
The additional Haldane pseudopotentials and the total 
Haldane pseudopotentials are listed in Table |VI| for the 
most common case of spin-1/2 fermions. For the case 
of spin-1/2 fermions, only odd m pseudopotentials are 
relevant for the Hilbert space where this effective Hamil- 
tonian acts. For a more general case of fo filled flavors, 
the pseudopotials V2(m, /o) are. 



V2(m,/o) = y2(m) + /„Vf(m), 



(32) 



where the coefficients V2{m) are given by Eq. (10) and 



VI 



are listed in Table [nj while Vf (m) are listed in Table 
It is interesting to note the significant difference between 
the corrections to the pseudopotentials for < < 1, 
listed in Table [III ^^'^ those for spin polarized states with 
1 < < 2, listed in Table VI For example, the correc- 
tion to the m = 1 pseudopotential, which is crucial in 
determining the gap of the Laughlin type states, is nega- 
tive and about four times larger in magnitude compared 
to the < I' < 1, indicating a higher reduction of this 
pseudopotential in the filling factor range 1 < v < 2. 

The three-body interactions between partially full fla- 
vor electrons remain unchanged and given by Eq. ( |11[ ), 
thus, in the spin-1/2 case the three -bo dy pseudopoten- 
tials would be those listed in Table HI for the S = 3/2 



states of Eq. ( A3 ), constructed for the partially filled fla- 
vor with spin ^. 
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TABLE V: CoefScients of the second order correction to the 5 = 3/2 and 5=1/2 three-body Haldane 
pseudopotentials in the n = 1 LL. For S = 1/2 and 4 < m < 6 the pairs hsted have been orthonormahzed by 
rotating only the state with a — 2, as described in Appendix [Al 



m 


3 


5 


6 


7 


8 


9 


V3/2{k'l',kl) 


(0,1) 
-0.0147 


(1,1) 
-0.0054 


(0,2) 
-0.0099 


(2,1) 
0.0005 


(1,2) 
-0.0009 


(0,3) (3,1) 

/ -0.0088 0.0007 \ 
I 0.0007 0.0033/ 


m 


1 


2 


3 


4 


5 


6 


{o,k,l) 
Vi/2(o''k'l' , akl) 


(1,0,0) 
-0.0319 


(2,0,0) 
-0.0305 


(1,1,0) 
-0.0131 


(1,0,1) (2,1,0) 

I -0.0009 -0.0004 \ 
I -0.0004 -0.0100 I 


(1,2,0) (2,0,1) 

( -5x10"^ -0.0056 \ 
V -0.0056 0.0229 / 


(1,1,1) (2,2,0) 

/• 0.0067 0.0017 \ 
1 0.0017 -0.0010.) 



TABLE VI: Coefficients of the second order corrections to the two-body Haldane pseudopotentials for spin- 1/2 
fermions with a state with completely full majority spins and partially empty minority spins in the n — Q LL. 



m 





1 


2 


3 


4 


5 


Vf(m) 
V2(m,l) 


-0.3662 
-0.7119 


-0.0959 
-0.1287 


-0.0268 
-0.0380 


-0.0078 
-0.0133 


-0.0023 
-0.0056 


-0.0007 
-0.0029 



V. SUMMARY AND DISCUSSION 

We have derived effective Hamiltonians which account 
for quantum fluctuations in full and empty Landau lev- 
els (Landau level mixing) to leading order in perturbation 
theory for the cases of a partially filled n = Landau level 
and a partially filled n — 1 Landau level. These effective 
Hamiltonians describe fractional quantum Hall physics 
in the filling factor ranges Q < v < 2 and 2 < v < A 
respectively in the case of spin-1/2 fermions. In both 
cases the effective Hamiltonians are a sum of two- and 
three-body terms. The three-body terms are responsi- 
ble for particle-hole symmetry breaking within the Lan- 
dau level. There has been considerable interest in these 
quantum fluctuation corrections because they are likely 
to play a decisive role in several outstanding problems in 
the fractional quantum Hall regime. 

Our work has been developed essentially in parallel 
with an analytic study by Peterson and Nayak,!^ and 
a numerical study by Simon and Rezayi.l^ The three 
works toghether provide a comprehensive view of the 
leading perturbative Landau level mixing corrections to 
the effective Hamiltonian in the lowest and second Lan- 
dau levels, and they complement earlier works by Rezayi 
and Haldane,!^ Murthy and Shankar,^ and Bishara and 
Nayak.l^ The source of discrepancy on the two-body 
pseudopotentials with the work of Bishara and NayalP^l 
is now understood. In Ref. '9', this discrepancy was shown 
to arise from the normal ordering that must be kept in 
the three body interactions. When this normal ordering 
is dealt with properly, our two- and three-body interac- 
tions are exactly the same as those derived in Ref. ,9, 

Our two-body pseudopotentials in the lowest and sec- 
ond excited Landau levels are in agreement with Refs. [9| 
and|43l In the lowest Landau level the two-body pseu- 
dopotentials are also in agreement with the results of an 
earlier study of the perturbative expansion of the spec- 



trum of hydrogenic atoms in a strong magnetic-field.'^ 
Our three-body pseudopotentials are in agreement with 
those of Refs. [9l and [42l and with those of Ref. |43]in the 
lowest Landau level. There is an appreciable difference 
with the three-body pseudopotential differences in the 
second Landau level reported in Ref. US} which is likely 
to arise from finite size effects in this numerical study. 

Now we would like to comment on the range of va- 
lidity of the perturbative approach. It is impossible to 
know the range of validity of any perturbative expansion 
without a sense for the relative size of the higher order 
corrections. We do not believe it is feasible to carry out 
the same type of analysis we have discussed in this paper 
to higher orders in n for the many-body problem in an 
analytic fashion. The reason is that it not possible to 
construct effective Hamiltonians to order or higher, 
projected onto the Landau level of interest, without ex- 
plicitly computing the energies to order k. In other words 
it is necessary to solve the many-body problem exactly to 
order k to be able to construct an effective Hamiltonian 
to order projected into the degenerate manifold. 

In spite of this seemingly insurmountable task, it is 
possible to get a sense, at least heuristically, for the size 
of the higher order corrections. One way this can be done 
is by studying the two-body problem to higher order in 
K. This has been done, indirectly, for the lowest Landau 
level in the context of the problem hydrogenic atoms in 
a strong magnetic field in Ref. HOI More specifically from 
the coefficients a''^ listed in Table II of Ref. i5Cil one 
obtains the following expression for the energies of two 
particles with relative angular momentum m in the lowest 
Landau level to order in units of Wc, 

= E"o!i (-^«)' + 0(«'). (33) 
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By reading the values of the coefficients ;„ from Ta- 
ble II in Ref. UHl one learns that a'^l^ decreases by about 
an order of magnitude at every ordei^^. This indicates 
that higher order corrections remain parametrically small 
even at k ~ 1. This observation suggests that higher 
order corrections might remain small even at values of 
K ^ 1 in the Lowest landau level, and specially so in the 
dilute limit where the energies of the two-body problem 
are expected to dominate. 



It is not possible to directly extract the energies of two 
fermions in the second landau level togheter with the 
completely filled lowest Landau level from the results of 
the bare two-body problem of Ref. [501 This is because the 
bare two-body problem fails to account for basic many- 
body effects like Pauli blocking that already arise at order 
K in the energies in the second Landau level. It is thus 
hard to anticipate at this point the relative size of the 
contributions in the second Landau level. 



Most fractional quantum Hall samples have similar 
electron densities. For this reason the external magnetic 
field strength tends to be smaller, and k correspondingly 
larger, for the experiments in the second Landau level 
than experiments in the lowest Landau level. Addition- 



ally by comparing Tables pi] and IV it is evident that even 
at a fixed field, the two-body pseudopotential corrections 
are larger for the second Landau level than for the lowest 
Landau level. Quantum fluctuations in Landau level oc- 
cupations are therefore more likely to be important in the 
second Landau level case in which the fractional quantum 
Hall effect can be enriched by the appearance of striped 
states and even-denominator incompressible states. 



Finally we would like to connect our study to the prob- 
lem of the nature of the incompressible state observed at 
filling fraction v — 5/2 in GaAs. For the spin polar- 
ized case, the Moore-Read Pfaffian is known to be the 
unique highest density zero- energy state of a repulsive 
three-body HamiltoniarP^I^SHIZl for which only the low- 
est angular momentum three-body state is energetically 
penalized, namely l3/2('Ti = 3) > and the remainder 
pseudopotentials vanish.!^ Conversely the anti-Pfafffian 
is expected to be the ground for the particle-hole conju- 
gated Hamiltonian,'^ which has 1^3/2 = 3) < 0. We 
have found in agreement with Refs. P and 3^ that the 
leading value of this pseudopotential is negative. Never- 
theless the two-particle pseudopotentials, which we have 
corrected for in this work, have been found to play de- 
cisive role in this competition, and they could drive the 
system into a compressible phase.'^SHUl A reliable assess- 
ment of the influence of LL mixing on fractional quantum 
Hall states at even denominator fractions thus awaits the 
application of our pseudopotentials in many-body exact 
diagonalization studies. 
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Appendix A: Three-body states and 
pseudopotentials 

In this appendix we outline some properties of the 
three body states and their associated generalized Hal- 
dane pseudopotentials. We begin with the n = LL. 
We are interested only in translationally invariant poly- 
nomials describing the internal state of relative motion of 
the three particles, and thus omit polynomial factors in 
the center of mass coordinates (zi + Z2 + z^)/^. Transla- 
tionally invariant polynomials would depend only on two 
independent translationally invariant coordinates, which 
can be chosen as 



The permutation operators on these coordinates act 
as a reflection, and as a composition of a reflection and 
rotations by ±27r/3,^ 



H12- 



1 

-1 




where the matrices are understood to act in a colum vec- 
tor of the form ( 1^ ) . 

The fully spin polarized S = 3/2 states, whose poly- 
nomial part is fully antisymmetric, were constructed in 
Ref. [511 and read as 



3/2 



1 



2\fe 



kl ~ rys,i-i (^a + ^b) 

^kl 



(Za + iZbr ~ {Za ~ iZb) 



2i 



(A3) 



with Z^f = 23'+2'=+i[7r3(3?-t-A:)!fc!]i/2, k>0,l>l. 
can be expanded as. 



/2 

kl 
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*fei 73/2 J- 



3/2 J m-j 

y3/2 jkl ^fc^O 



(A4) 



where m = 2fc + 3Z is the orbital angular momentum of 



, and 



m m—2 



= E E - J - / , J>12|/>3, (A6) 

j=0 j'=0 



where |m, »Ti')ij abbreviates for the state with center of 
mass angular momentum m and relative angular momen- 
tum m' for particles i and j in the n = LL, and \l)i for 
the state of particle i with angular momentum I in the 
n = LL. Cki.jj' can be found to be: 



A: 31 

Jkl "EE 

p=0 g=0 



sin(q7r/2)(5j,5+2p- (A5) 



In order to evaluate the pseudopotentials from Eq. ( 12 ) 
it is convenient to decompose the state 5"'^;^ into products 
of states with well defined relative numbers for the pair 
of particles 1 and 2, and states with well defined numbers 
for particle 3, as follows: 



C, 



kl,jj' 



^ki \ ■> 



(m-j-/)!j!j'!. (A7) 



Consider now two states, and with the same 

angular momentum, m = 2fc + 3 = 2k' + 3^'. Their 
associated generalized Haldane pseudopotential matrix 
elements, computed from Eqs. (12) and (15), is 



v,/^{k'i',ki) = 6j2J2 ^fc'''."' E E ^i^w E 

j=0 j'=0 n=l 



1=0 i'=0 



.(0, i\vi2\n, i + n)rr{n,j + n\vi2\0,j)r j^m+n-j' , . o^ 

2«„ -«j+nj' ^i+n,i' ' \^^) 



where i?^ „j, are given by Eq. ( C4 ) , and the Coulomb 



matrix elements by Eq. (C6). Equation (A21) can be 
used to obtain the S = 3/2 pseudopotentials listed in 
Table Uni 

Let us now discuss the S = 1/2 three-body states. 
We begin by considering the bosonic analogue of 
namely 



kl 



{Za + iZi,y' + (Za - iZb) 



311 



, (A9) 



with k > 0, I > 0, and we have not normalized these 
states yet. These states form a complete orthogonal basis 
for the fully symmetric translationally invariant polyno- 
mials. They can be expanded as 



polynomials with a mixed permutation symmetry, repre- 
sented by the following Young tableaiP^J 



1 


3 


2_ 





These primitive polynomials ar^^ 

/3l = Zb, 132 = ZaZb- 



(A12) 



(A13) 



The most general polynomial with the symmetry of the 
tableau ( A12| is a product of either one of these polyno- 
mials with the most general symmetric translationally 
invariant polynomial $fci, namely. 



with 



2k+3l 

<fk, = V f^Q ztz 



2k+3l-j 



(AlO) 



k 31 

Jkl "EE 

p=0 q=0 



cos(q7r/2)(5j,5+2p. (All) 



To construct the spatial part of the S' = 1/2 states, we 
start by writing the lowest order translationally invariant 



with a 



2\k 



yl/2 

^akl 



{zi + zt) 



{Za + iZbf + {Za - iZbf 



(A14) 

{1,2}, fc > 0, ^ > 0, with the normalization 



constants Z\'j^i given in the text below Eq. (|18|). This ba- 



sis is complete but not orthogonal. Fortunately, the inner 
products can be computed analytically and are listed in 
Eq. ( 20 ) , thus othonormalization becomes trivial. 

The spin part of the S* = 1/2 wavefunction h as th e 
symmetry of the Young tableaux conjugate to (|A12| ). 
Therefore, the fully antisymmetric wavefunction ia 



12 



(A15) 



where A — (1 — n23+ni2n23)/\/3 is the antisymmetriza- 
tion operator between f and | particles, and Ilij is the 
unitary operator corresponding to elementary permuta- 
tions between i and j. 

The three-body interaction in Eq. (Ill is rotationally 
invariant, hence, its associated three-body pseudopoten- 
tials are diagonal in the relative orbital angular momen- 
tum of I'^l^Zi), m = a + 2k + 31, therefore, for any two 
states with the same angular momentum we can write 



where m = a + 2k + 31, and the indices {/ii^7} stand for 
any permutation of {123}. The decomposition is gener- 
ally dependent on the particle ordering since the orbital 
part alone of the S' = 1/2 states, '^1^^, is not an eigen- 
state of all permutations. C^^; ^j, can be found to be 



^<7kl,jj' 



m — J 
f 



i-iy 



2m+j' 



{m-j-f)ljlfl, (A19) 



Vi/2{a'k'l',akl) 



n 



^l^ki) 



with 



2%',U'{ho + hl + h^)^^J^^, (A16) 



where n = hi + n2 + fis, V±_ is the projector into the 
othogonal complement to the n = eigenspace. In the 
second line of Eq. ( A16 1 it is understood that only the or- 



bital part of the wavefunctions is involved, and it follows 
from Eq. (A15| using the fact that interactions are spin 



independent. The computation of Haldane pseudopoten- 
tials reduces to the computation of the matrix elements 
of /iQ and hi, which stand for 



ho ^ -vi:iVA_\'P^vi2, hi = -vi3Vj_^Va_V23- (A17) 
n n 



In analogy with Eq. (A6), we can decompose '^J/^i as 



akl 



EE^'^'u.'I"^-J-^"'JVI/)7- (A18) 



9crjkl 


„ fl/2 

~ Jj~l,kV 






9ijki 


2 






9ijki 


f 1/2 






92jkl 






Jj,kU 




4 




92jkl 






Jj.kl) 




4 





(A20) 



These expressions need not be derived independently, but 
can be obtained by deriving only the decomposition cor- 
responding to {/ii'} = {12}, and then applying suitable 
permutation operators as repres ented in Eq. (A2). With 
this the evaluation of Eq. (A16l leads to. 



m m—t 



m m—j 



^{0,i\vi2\n,i + n)rr{n,j +n\Vi2\0,j)r „m+n^i' T>m+n~j' 

i+n,j' j+n,i' ' 



i=0 i'=0 
■m rn — i 



j=0 ]'=0 n=l 
m m—j c 



2"n 



^'?.,i,hi^'{' 

a' k' L' J- akl 



^ ^ (-l)'Cl3j,,,, , ^ ^ C^l^ ., ^ r(0, ^1^121", ^ + n)rr{n, j + n\vi2\Q,i)r ^„,+n^,' ^m+n-j' 



j=0 i'=0 



j=0 j'=0 



n=l 



2"n 



i-\-ri,j' j-\-ri,i' 



(A21) 



Notice that all the matrix elements are purely real. Com- 



bining these expressions with Eqs. ( C6|), (|C4| and (|A16|), 
one obtains the values listed in Table [ 



II 



The derivation in the n — 1 LL goes through in a com- 
pletely analogous fashion. By raising the kinetic energy 
of the three particles, with the operator a\a}2a\, we easily 
obtain the required representation of the states in one- 



to-one correspondence with the n = LL. 

To obtain the S* = 3/2 three-body pseudopotentials in 
the n = 1 LL we apply a\a\a\ to Eq. (A6). Using the 
first quantized version of the three-body interactions ap- 



pearing in Eq. (27), we arrive at the following expression 



for the pseudopotentials in the n = l LL, 



13 



m rn—i 



j=0 j'=0 '-n=l 



i=0 i'=0 



r{0,i\vi2\n, i + n)r,.(0, J>i2|n, j + n)r ,n+n-i' j^m+n-j' 
2nji i+n.j' ^j+n,i' 



r{2,i\vi2\l,i- l)rr{l,j - l|fl2|2, j),. pm-1-/ 



(A22) 



where we have introduced the notation, 

r(0, i\vi2\n, i + n)r = ^(0, i|wi2|n, i + n)rR'^'ji^ ■ ■ ■ 

r{2, i\vi2\n + 2,i + n)rR';+li- (A23) 



From Eq. (A22) one obtains the pseudopotentials for 
5 = 3/2 states in the n = 1 LL hstcd in Table |v) 



To obtain the 5 = 1/2 states in the n — 1 LL we apply 



a\ala\ to Eq. (A18). From the interaction appearing in 



Eq. (271 we obtain the analogue of Eq. (A16I, now with 



ho and hi replaced by 



ho = -Vi3V± , ^ „ 'P±Vi2, hi = -Vi3V±^r^—V±V23- 

n — 6 n — 6 



(A24) 



The matrix elements for ho are, 



m ra—i 



'^a'l'l'^o'^akl ^ ^ E E ^^^k'l',ti' E E ^^kl.jj' E 

j=0 ]'=0 '-n=l 



i=0 i'=0 



..{0,i\vi2\n,i + n)rr{0,j\vi2\n,j + n)r rn+n-i' j^m+n-f 

i+n,j' ^j+n,i' 



r{2,i\vi2\lj - l)rr(l,i - l|fl2|2, j),. ,n-l-i' pm-1-/' 



(A25) 



The expression for 
ter replacing C^^ 



a' k' L' 



(-132 



is the 
and 



same af- 



(-l)'C^?fc,;,^,,, in analogy with Eq. ( |A21[ ) for the n = 
LL. The 5 = 1/2 Haldane pseudopotentials in the n = 
1 LL, listed in Table |Vj are given Vi/2{(j' k' I' , a kl) = 



Eq. (|A16|) for the n = LL. 



(h, 



which is the analogue of 



In spite of how cumbersome Eqs. (A22| and (A25) 



might look, its evaluation is very efhcient, and each pseu- 
dopotential takes only a few seconds to evaluate using 
Mathematica in a conventional laptop computer. 

We listed in Tables III and [V] pseudopotentials matri- 
ces up to 2 X 2 in size. For this case the two-dimensional 
subspaces for given m contain one state for a = 1 and 
one for (7 = 2. We have orthogonalized these states by 
rotating the a = 2 state only. Let us call the matrix 
in the non-orthogonal basis V, which is obtained from 
Eqs. ( |A16[ ) and |A21[ The matrix listed in Tables [m 
and [y] corresponds to (B-^)^^VB^^, with B the change 
of basis matrix, 



B 



1 (*l|4'2) 

(*l|*2)2 



(A26) 



{^i\'^2) is a shorthand for the overlap between cr = {1, 2} 
states, appearing in Eq. (20). 



Appendix B: Effective interactions in the n — 1 LL 

In this appendix we describe in more detail how the 
interactions in the n = 1 LL were obtained. Consider 
the effective Hamiltonian to order in the n = 1 LL, 



"^2 = "Hi - ^ ^ Wl2,34t'56,78 X ' ' ' 

PoclAf. C7,.C8uP± j^^ P-i- 4 A 4<T C3<T C4 A Po , (B 1 ) 

where the integers abbreviate single-particle kinetic and 
guiding center quantum numbers (e.g. {1 <^=> nimi}), A^i 
is the number of particles in the partially filled n — 1 
LL, Po is the projector into the many-body eigenspace 
with h = Ni, and P± the projector into its orthogonal 
complement. 



The second order term in Eq. (Bl I can be viewed as a 



sequence of scattering process in which a pair of particles 
is taken from states "3" and "4" , contained either in the 
completely filled n — LL or the partially filled n = 1 
LL, and placed into states "1" and "2" with higher total 
kinetic energy. Subsequently, the particles are removed 
from these virtually excited states by operators "7" and 
"8", to be finally placed back into states "5" and "6", 
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which are contained either in the n 
LL. 



LL or the n = I 



Threrefore any term with non-vanishing matrix ele- 
ments in the h = Ni eigenspace must have the outer- 
most destruction and creation operators (i.e. Cg^^, Cg^, 
csa- and c^x) with Landau level indices either or 1. 
Therefore, the kinetic energy of the "incoming" parti- 
cles, i.e. 713 + "-4, or the kinetic energy of the "outgo- 
ing" particles, i.e. + n-ej is allowed to take only the 
values {0, 1, 2}. It can be shown that there is a kinetic 
energy balance condition between "incoming" and "out- 
going" labels, namely the terms with ns + ^ + uq 
vanish, and they correspond to scattering type diagrams 
that have tad-pole and self-energy insertions in their 
legs. In summary, we can label the allowed interac- 
tions by the total incoming/outgoing pair kinetic energy 
"3 + = "5 -I- 716 = {0, 1, 2}. 

As an example we discuss in detail how to obtain the 
interaction corresponding to incoming pairs with total 
kinetic energy 713 -|-n4 = 2, labeled V^^(2), and will leave 
the verification of the remainding terms to the interested 
reader. This term corresponds to all the two-body terms 
arising from Eq. (Bl) with 71,3 = 714 = 715 = rig = 1, 
after normal ordering of the operators is performed. In 
order to have a non-vanising contribution the creation 



operators and must raise the kinetic energy of 
the incoming pair, thus 711 + ri2 > 2. Similarly we have 
that 777-1-778 > 2. Since the 77 = LL is assumed to be 
completely full, these labels are additionally constrained 
to satisfy 77,; > 1, for i = {1, 2, 7, 8}. 

With these constraints we can see that there are two 
possibilities for these middle operators. Either neither 
of them corresponds to the partially filled n = \ LL, i.e. 
rii > 1, or else only one creation and only one destruction 
do correspond to n = 1 LL, for example 772 = 777 = 1 and 
ni = 778 > 1- For the first possibility we see that the 
indices {1,2,7,8} must be fully contracted in order for 
the virtually excited pair to go back into the 77 = 1 LL. 
Let us name the terms arising from this first possibility 
as Va^{2). We can write 



^2 

(2) = -— '"l2.,34V5ejs{c7f,Cs^c\^cl„) X ••■ 



1---8 



;^^:p^^C5,Cg,c3.C4., (B2) 



Vb(2) 



K 



Y Wl2,34f56,78cL4M 



1---8 



-771 — 1 

[{Csuc\x)c7^,cl„ - (C7,,ci^)c8^4^)+{1 ^ 2} 



C3ctC4A, 

(B3) 



where again we understand the Landau level index of the 
uncontracted operators to be in the n = 1 LL. 

It is clear from Eq. (B3), that V;,(2) has a contri- 



bution to the effective three body interactions, but it 
also contributes to the two body interactions after due 
normal ordering of the operators is performed. The 

contributes to the terms in 
le other hand, the two body 



three body piece in Eq. ( |B3 
Eq. (l26|, with 777 > 2. On t 



term, which we label VP(2), makes up the remainder 



of the interaction V^''(2) appearing in Eq. (24 1, namely 
V2^(2) = Vf (2) 



(2). 



A very similar analysis leads to the forms for the 
remainding two- and three-body interactions given in 



Eqs. (241, (251, (26), and (31). 



Appendix C: Two-body states and Haldane 
pseudopotentials in the 77 = 1 LL 

In this section we discuss several useful properties of 
the two-body states in the n = 1 LL and their associated 
Haldane pseudopotentials. 

We begin by discussing two-particle states. Consider a 
state for particles i and j in which each of them has well 
defined kinetic and guiding center quantum numbers: 



|r7i7r7i)i|7727772) 



' ' ^ |n)^ 

%/77l!772!77il 17772! 



(CI) 



where and are th e kinetic energy and guiding center 
raising operators .l^^^^ Although we often loosely refer to 
the guiding center quantum number 777^ as the angular 
momentum of particle 7, notice that the physical angular 
momentum is actually 777^ — 77^ . 

An alternative representation for two-particle states is 
obtained by constructing states with well defined center 
of mass and relative coordinate quantum numbers, for 
both the kinetic energy and guiding center labels, as fol- 
lows 



where now we understand the Landau level index of the 
uncontracted operators to be in the n = 1 LL, and ^{n) 
restricts n > 2. 

For the second possibility an extra pair of uncontracted 
creation and destruction operators in Eq.(Bl) will have 
indices in the 77 = 1 LL. Let us name this term Vf,(2), 
which reads as 



I iVM)c 1 77777) r = 



VN\Mln\m\ 



|0), (C2) 



where the four operators A = {ai + aj)/V2, B = (6^-1- 
hj ) I \/2 , a = [ai — aj)/ \f2 and h = {bi~-bj)/ \/2 , commute 
with each other. These states are related to those in 
Eq. (iCll) by the unitary transformation 
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ni+n2 rni+m2 

i/=o fj.=a 

\ni+ na - v,mi + 1712 - li)c\v^i)r, (C3) 

where R^-^ is a real, orthogonal, and symmetric L x L 
matrix. The explicit form of i?^^ is 



min{m' .m) 

E 

/i — 7\/aa;(0,m+m' — L) 



L ~ m 
m' — ^ J \fi 



(-l)^ (C4) 



The properties of imply that the inverse transfor- 



mation is still given by Eq. (CI) with the exchange of 
roles {i c,j r}. 

For two isolated particles in vacuum the center of mass 
and relative basis of Eq. ( C2 1 is the most natural basis 



to find eigenstates of the two-body problem in the pres- 
ence of strong magnetic fields. The reason is that the 
interaction is almost diagonal in this basis: 



,{N'M'\r{n'm'\v,j\NM)c\nm)r 

= r{n'm!\vij\n'm)r5rn-n,rn'-n'5NM'5MM' ■ (C5) 



For the Coulomb potential, Vij — l/|rj — rj|, with 
lengths measured in 1^ units, the non-vanishing matrix 
elements are: 



Equation (C6) can be obtained by combining the explicit 
expressions for the relative coordinate wavefunctionJ^Sl 
with the useful integrals computed in the Appendix of 
Ref.lMl 

For the iV = /A^ll + 2-body problem, in which the 
lowest Landau level is completely filled and only two par- 
ticles are in the n = 1 LL, the basis of pure relative and 
center of mass quantum numbers is no longer the most 
convenient choice to describe the relative motion of these 
two particles. The reason is that generally states with 
well defined relative kinetic energy, n > 1, will have a 
probability amplitude for the individual particles to oc- 
cupy the n = LL, which is forbidden by the Pauli ex- 
clusion principle. The most natural states would look 



entangled in both the basis of Eq. (CI I and Eq. (C2|, 
and have the form 



|M,m)i = al4|M,0)c|m,0) 
1 



V2 



(|M, 2)e|m, 0), - \M, 0)c|m, 2),). (C7) 



These are the states that d eterm ine the Haldane pseu- 
dopotentials in the n = 1 LlP^EH appearing in Eq. ( 29 1 



To simplify notation, we have denoted them simply in 



Eq. (29) by |m)i, since the pseudopotentials are diago- 
nal on and independent of the center of mass label M . 

The computation of Haldane pseudopotentials in the 
n = 1 LL is simplified by transforming the two body in- 



teractions of Eqs. ([24|), ([25|), (|26|, and (|31|, into their first 
quantization version, wE 



rich we describe in the remainder 



of this Appendix. For V'^^{2) appearing in Eq. ( |24[ ), we 
have the equivalent first quantization form 



I ji + m 

tt{1' + \j\)U'W. 



3F2 



2|J|! 

-MjI 



1 1 

|j| + l,l/2~V'^ 



(C6) 



where m > 0, m' > 0, n > 0, n' > 0, m' — rt' = m — n, 
J =m-n,l = n+{j- |j|)/2, and 1' = n' + (.? - |.?|)/2. 



(2) = -n- 



El 1^2,) (1,2, 



1,2 



0{ni,n2) 
rji + — 2 



(C8) 



where |li2j) is a shorthand for the state \nimi) i\n2'm2) j , 
and 9(ni,n2) restricts rii -|- 77.2 > 3 and ni > l,n2 > 1. 
Combining Eqs. and ((GTI) with Eq. (IC8| leads 



to the following semi-analytic expression for the Haldane 
pseudopotentials associated with V^''(2): 



V 



2b 



1 °° 1 

(2)M^-iHVf(2)|m)i = -E2n 



(n, m -I- n|-(;i2|0, m)^ + r{n + 2, m + n|-yi2|2, m)1 



(^r{n, m + 7i|i;i2|0, m)r\/ {n + 2){n + l)/2 - r{n + 2,m + n|wi2|2, m)^ 



(C9) 



For V^^(O), appearing in Eq. (24), we have the equiv- 
alent first quantization form 



■ 1.2 



|li2j) (l,i2j |(5ni,0^n2,0 



v,j. (CIO) 
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This form leads to the following analytic expression for 
the Haldane pseudopotentials: 



V'\2){m)^-^,{m\Vt^{2)\m),= 

- ^,(2,m|i;i2|0,m-2)2. (Cll) 

We convert the remaining two-body interactions into 
first quantization using a slightly different technique, 
along similar lines to the original computation of the con- 
ventional Haldane pseudopotentials in the n — 1 LLpSl 
We write the interaction as a Fourier sum 



(C12) 



where abbreviates iMX^q- 

For Vf''(l) the interaction matrix appearing in Eq. ( 25 ) 



can be seen to correspond to the following first quantized 
interaction: 



J2 ^^'''''''' (6|e"?^-'-|5)(5|e-'^^-n6) 



5,6 



(C13) 



where |5) abbreviates single-particle state In^m^) (simi- 
larly for |6)), r(n5,ng) restricts ng > 1 and = 0, and 
/ is the spin multiplicity. Then using the properties of 
the matr ix elements of the single-particle density opera- 
tor™ 

one finds 



2 \^Qn{q)? ig.(r 

n=l 



(C14) 



It is interesting to note that this expression is nothing 
but the term of the well-known RPA screened poten- 
tial for a completely filled n = LL. The static screened 
RPA potential is, 



(C16) 



with Xq the static density- density response function. For 
a completely filled n = LL, the static density-density 
response function i^^ 



X„ = 



f 



fn fn 



^ n 



In ^ — ' n 



27r ^-^ n — n' 

(C17) 

where the Fermi occupation factor, /„, is 1 for n = and 
otherwise. The density form factors J^nn' (q) are given 
in Eq. ( |C21[ ) below. It is clear from Eq. ( |C17[ ), that the 
term of v^^^"^ is equivalent to Va\ji^) appearing in 



Eq. (CM) 



V^^l), appearing in Eq. (251, can be seen to corre- 



spond to the following first quantized interaction. 



91,92 
-iqi-rj 



g»9i-n X 



E |5.)(5. 



6,n6>l 

g-»92-r, + {j ^ jj^ (C18) 



where |6i)(6i| operates only on particle i, and |5j)(5j| 
only on j. It is not transparent in this representation 
that V^^j (l) has translational and rotational invariance. 
These properties can be made manifest if we write the 
position operators, r, in terms of the mechanical mo- 
mentum, TT, and guiding center, c, coordinates, from the 
following equations: 



where iJ'oniq)]'^ = |q|2"e"l«l /V(2""-!) is the modulus 
squared of the d ensity form factors between the an d n > 



Landau level^^^l^ (explicitly presented in Eq. (C21) 
below), and = / (Pq/{2n)-^. The Haldane pseudopo- 
tentials are then found to be 



Vf (1)(™) I dq qvl[L,{qy2)rLUq') x 



TT : 
C 



P 

r - 



'A(r), 
c 

Z X TT. 



(C19) 



The mechanical momentum and guiding center coor- 
dinates are related to the lowering operators as a = 
(tt, + iny)/V2, b = (c, - ic^y2.-48 49. 

The interaction in Eq. ( |C18[ ), projected into the n = 1 
LL, can be written as 



e-3«'/2 ^ ^^(^2/2) _ Ln{q'/2)] , (C15) 

where 7 « 0.5772 is the Euler-Mascheroni constant, 
Ei{x) = — dte~*/t is the exponential integral, L„ 
are the Laguerre polynomials, and Ln is the natural log- 
arithm. We substituted the explicit form of the Coulomb 
potential, Vq — 27r/g, in this expression to obtain the 
numbers listed in Table ITVl 



V,%(1) 



j(9l-92)-(ci-Cj) 



91 92 



J^ioi-qi)J^oiiq2) 



E 



-^ln(-g2)-^nl(gl) 

n 



(C20) 
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where J^n'nil) = (n'|exp(ig • z x 7r)|n) are the density 
form factors, given bj^^iMi 




V /t : \ V z / \ ^ 

(C21) 

for n' > n, and J-n'n{—<l) = J'mi'il)- Combinig this with 
the analogue relation for the guiding center coordinates, 



,(0,m|e^«-(^--=^)|0,m), = e-^L„,(|<7n, (C22) 
one finds the Haldane pseudopotcntials for V^'' listed in 



Table HVl 

For Vc^(l), appearing in Eq. (25), following an analo- 



gous procedure as the one just outhncd for V^^(l), one 
finds the following first quantized representation for the 
interaction projected into the n = 1 LL: 



91 92 



e-9^x«iJ-oi(-g2)-Fu(-gi] 



E 

n=l 



(C23) 



Combining Eq. (C23) with Eq. (C22) one obtains the 
Haldane pseudopotcntials listed in Table IV for V^^{1). 



Finally, the pseudopotcntials corresponding to the ef- 
fective interaction in the spin polarized n = LL, ap- 



pearing in Eq. (31 1, can be obtained similarly to those 
for V^''(l), and read as, 



fo 



2n 



e-39V2 + Eiiqy2) - Lniqy2)] , (C24) 



where fo is the number of fully occupied flavors {fo < /)■ 
The numbers associated with this interaction are listed 
in Table EH 



As an independent consistency check, we have com- 
puted the Haldane pseudopotcntials for the n = 1 LL 
appearing in Table |IV[ for the leading spin triplet states 
(odd m), in a more direct numerical approach using the 
two-body interaction matrices that appear in Eqs. (24) 



and ( 25 ) , and verified they coincide with those computed 



through the semi-analytical formulae discussed in this 
Appendix. 

We have presented our formulae in a manner that can 
be easily modified to incorporate cases other than pure 
Coulomb interactions, with the only requirement of ro- 



tational invariance. Only Eq. (C6) makes explicit use of 



the form of the Coulomb interaction. 
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